Improved setup for producing slow beams of cold molecules using a rotating nozzle 
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Intense beams of cold and slow molecules are produced by supersonic expansion out of a rapidly 
rotating nozzle, as first demonstrated by Gupta and Herschbach 1, 2]. An improved setup is 
presented that allows to accelerate or decelerate cold atomic and molecular beams by up to 500 m/s. 
Technical improvements are discussed and beam parameters are characterized by detailed analysis of 
time of flight density distributions. The possibility of combining this beam source with electrostatic 
fields for guiding polar molecules is demonstrated. 

PACS numbers: Valid PACS appear here 



I. INTRODUCTION 

Cold molecules currently attract great attention due 
to their potential application for precision measure- 

1 IfL 



ments [H, |3| , quantum information schemes [3, [Bj , de- 
generate quantum gases with complex interactions J, @] , 
and for cold chemistry 0, 0, H B E3, EI- Besides the 
approach of forming cold molecules out of a sample of 
ultracold trapped atoms [1, [H, 03, recent exper- 
imental efforts mainly focus on developing techniques 
of slowing or filtering, cooling, and trapping molecules 
produced in effusive or nozzle beams. Well established 
techniques are buffer gas cooling 0, [H[ , the deceleration 
of beams of pol ar molecules using time-varying electric 
fields 0, [H, , as well as filtering of polar molecules 
out of an effusive source using static or time- varying elec- 
tric fields [HI, [Til . Recently, the method of decelerating 
supersonic molecular beams has been extended to mag- 
netic [2fl |2~0 ] and optical fields [22j . Alternative routes 
to producing cold molecules have been demonstrated uti- 
lizing the kinematics in elastically or reactively colliding 
molecular beams fiol. I2I I2H] . 

Furthermore, producing beams of slow and cold atoms 
and molecules by mechanical means has been demon- 
strated P, 0, Hi| . In particular the technique of trans- 
lating a supersonic jet to low longitudinal velocities by 
means of a rapidly counter-rotating nozzle, pioneered by 
Gupta and Herschbach, holds the promise of producing 
cold, slow and intense beams of nearly any molecular 
species available as a gas at ambient temperatures [l], [|J . 
Attempts to improve the original arrangement have been 
made by the groups of H.-J. Loesch in Bielefeld and of 
M. DeKieviet in Heidelberg [2||. By controlling the sense 
and speed of rotation of the nozzle as well as by using the 
seeded beams technique, tunable beam velocities rang- 
ing from thousands down to a few tens of meters per 
second are feasible. Typically, the translational and in- 
ternal temperatures in the moving frame of the slowed 
beam are in the range of a few Kelvin and a few tens of 
Kelvin, respectively. Due to the extremely simple con- 
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cept, this technique can be realized by a relatively simple 
mechanical device. Therefore this versatile source of cold 
molecules may be particularly useful for novel studies of 
molecular reaction dynamics at low collision energies in 
the gas phase or with surfaces @, H, H3, HI] . 

The principle of operation relies on the Galilei trans- 
formation of particle velocities in the moving frame of 
the rotating nozzle into the resting laboratory frame. As 
the nozzle moves backwards with respect to the molec- 
ular beam motion the beam velocity v s is reduced to 
vo = v s — v ro t- In this way, Gupta and Herschbach suc- 
ceeded in slowing beams of krypton (Kr) and xenon down 
to 42 and 59 m/s, respectively, at beam temperatures 
down to 1.5 and 2.6 K, respectively Q. In their proto- 
type setup they used a tapered aluminium barrel with a 
hollow bore as a rotor. The special aluminum alloy was 
chosen for its favorable ratio of tensile strength to mass 
density, allowing high rotational speeds. However, ma- 
chining such a rotor barrel requires considerable efforts. 
The gas injection was realized by inserting a polymer 
needle into the rotor. However, erosion of this type of 
sealing and the resulting gas leakage was found to be a 
serious limitation. The nozzle orifice was fabricated by 
gluing a stainless steel disk with a 0.1 mm aperture close 
to the tip of the rotor. In this case, swatting of slowly 
emerging molecules at high rotational speeds by the ro- 
tor arm as it comes back after one full revolution may 
be an issue. Furthermore, a general drawback of produc- 
ing slow beams using a moving source is the fact that 
the flux of molecules Nq on the beam axis drops down 
sharply with decreasing laboratory beam velocity Do and 
with increasing distance d of the detector from the nozzle, 
F(vo,d) tx (vo/d) 2 . This is a consequence of the beam 
spreading in transverse directions which becomes more 
and more important as vq decreases and hence the flight 
time goes up. The aim of the improved setup presented in 
this paper is to remedy the drawbacks mentioned above. 



II. EXPERIMENTAL SETUP 

The experimental arrangement used for characterizing 
the new rotating nozzle setup is schematically depicted 
in Fig. [T] The upper part represents the moving molecu- 
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FIG. 1: Schematic overview over the experimental arrange- 
ment used for characterizing the rotating nozzle molecular 
beam source. 



lar beam source including a hollow rotor and a skimmer 
with 1 mm diameter. The nozzle chamber is pumped by 
an 8000 1/s oil diffusion pump backed by roots and rotary 
vane pumps. Typical background gas pressure during op- 
eration with argon (Ar) at po = 0.3 bar nozzle pressure 
is pb g w 10~ 4 mbar. The molecular beam follows a tra- 
jectory through the following elements: At a distance 
of about 20 mm behind the nozzle it enters an interme- 
diate chamber through a skimmer with 1 mm opening 
diameter. The intermediate chamber is pumped by a 
600 1/s turbo pump providing a background gas pressure 
Pdetect ~ 5 x 10 -9 mbar during operation. The intermedi- 
ate chamber is connected to a detector chamber pumped 
by a 150 1/s turbo pump to maintain a background pres- 
sure pdetect < 10~ 8 mbar. 

For the purpose of guiding beams of polar molecules, a 
set of four parallel steel rods is installed around the beam 
axis between the skimmer and a commercial quadrupole 
mass spectrometer (QMS) (Pfeiffer QMS200) which we 
use as a detector (see Sec. IVT)) . The mounting discs used 
for holding the rods also serve as apertures for differential 
pumping between the intermediate chamber and the de- 
tector chamber. The total fight distance from the nozzle 
to the detector amounts to £ = 35 cm. Since the positions 



of the quadrupole guides and the QMS detector cannot 
be adjusted, the beam axis is initially aligned by sight 
by adjusting the position of the skimmer. For additional 
characterization of the beam, a Pitot tube is installed 
about 5 mm behind the ionizer of the QMS detector. It 
consists of a cylindrical chamber 4 cm in diameter and 
10 cm in length having a conical aperture on the side 
close to the QMS and a cold cathode ion gauge on the 
opposite side as an end cap. When a continuous molec- 
ular beam (nozzle at rest) is directed into the detector 
chamber, the residual gas pressure inside the Pitot tube 
rises until equilibrium between incoming flux and out- 
going flux of molecules effusing out of the Pitot tube 
aperture according to its molecular flow conductance is 
reached. Thus, from the difference between pressures in- 
side the Pitot tube and in the QMS chamber, Ap, the 
flux of molecules on the beam axis, TVo, can be inferred 
according to No = ApAav / '(4fceT) . Here A, a, and v de- 
note the area of the entrance channel into the Pitot tube, 
a scaling parameter to account for the channel length, 
and the mean velocity of the molecules at ambient tem- 
perature T . The on-axis beam density, which is obtained 
from the molecule flux using n = Nq/^vqA), is compared 
to the one determined by the QMS signals in Sec. IIVI 

The experimental realization of our rotating source of 
cold atoms and molecules is depicted in Fig. O(a). The 
following main modifications with respect to the work of 
Gupta and Herschbach have been implemented: The ro- 
tor is made of a thin carbon fiber (CFK) tube with 4 
mm outer diameter and 38 cm in length. CFK is supe- 
rior to any metallic material regarding the ratio of tensile 
strength to density by at least a factor of 4. Moreover, 
the small mass of the CFK tube makes it simpler to 
suppress deterioral mechanical vibrations that may oc- 
cur due to slight imbalance of the rotor at high speeds 
of rotation. The CFK tube is sealed on both ends by 
gluing thin-walled (0.25 mm) titanium grade 5 ferrules 
as end-caps into the inner surface of the CFK tube using 
special epoxy for CFK (Fig. [U (c)). One of the ferrules 
has a borehole d — 0.1 mm in diameter at a distance of 
0.5 mm from the closed end, thus forming a well defined 
nozzle aperture. A simple alternative that we also tried 
is capping the CFK tube by gluing circular pieces of cap- 
ton foil to front and back sides of the tube. Capping the 
CFK tube this way has proven to be stable enough to 
withstand the high centrifugal forces acting at rotational 
speeds up to f rot — 350 Hz (21000 rpm). The nozzle ori- 
fice was fabricated by inserting a 50 /jm wire half way in 
between the CFK front side and the capton foil prior to 
gluing. After the epoxy has hardened the wire is pulled 
out and an approximately 1mm long channel remains. 
This method allows to easily realize small nozzle chan- 
nels of diameters < 100 /im, which limits the gas load 
and thus the required pumping speed. Moreover, using 
this method of fabrication the distance of the nozzle from 
the tip of the rotor is kept minimum, thereby suppress- 
ing swatting of slow molecules by the rotor. However, 
the expansion out of a long nozzle channel is unfavorable 
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FIG. 2: Overview (a) and close-up views (b,c) of the exper- 
imental setup, (a) The nozzle orifice is mounted on one tip 
of the hollow rotor and produces molecule pulses along the 
beam axis each time it passes nearby the skimmer. The gas 
is injected into the rotor through a rotary feedthrough from 
below. The motor is mounted to a massive brass mount on 
top of the rotor, to which it is connected by a flexible coupler 
(silicone tube), (b) Sectional view of the central part of the 
setup, (c) Close-up view of the rotor tip, which contains a 
bored thin-walled titanium ferrule as a nozzle. 



compared to a thin- walled orifice as far as adiabatic cool- 
ing conditions are concerned. Besides, nozzles fabricated 
using this procedure vary considerably with respect to 
length and shape of the nozzle channels. For all measure- 
ments presented in this work, bored thin-walled titanium 
ferrules are used. 

The CFK tube is inserted and glued into the core of the 
rotor consisting of an aluminium or brass cylinder (Fig. [5] 
(b)). The latter has a borehole along its axis for insert- 
ing the bored shaft of a rotary feedthrough from below 
and the shaft of a flexible coupler attached to a motor 
(faulhaber-group.com) from above. In order to avoid any 
imbalance it is important to perfectly align the motor 



axis and the feedthrough axis. This is ensured by in- 
serting the two axes into the same through-hole through 
the rotor core. The rotor including the metal core is 
first balanced statically by suspending the rotor core be- 
tween two needles in the center of gravity and by abrad- 
ing CFK material from the heavier side of the rotor until 
it is perfectly balanced. In a second step, the rotor is 
dynamically balanced by adding screws to the rotor core 
at different places and by monitoring the vibration spec- 
trum using acceleration sensors attached to the setup. 
Using this procedure, vibrations at high rotation speeds 
resulting from the rotor are small compared to vibrations 
induced by the motor itself. Therefore, additional masses 
are attached to the motor (Fig. [U (b)), which reduce vi- 
brations at frequencies > 100 Hz. Gas injection into the 
rotor is realized using a commercial rotary feedthrough 
sealed by ferrofluidic metal sealings. The bored shaft 
of the feedthrough is inserted into the rotor core and 
O-ring sealed. Metal-sealed feedthroughs (ferrotec.com) 
features very high leak tightness (< 10~ n mbarl/s) at 
low friction torque. 

The rotor axis is aligned with the motor axis by con- 
necting the rotary feedthrough and the motor to a U- 
shaped brass construction. This brass mount is broken 
in two parts to introduce a rubber layer in between the 
two parts as a damping element for suppressing vibra- 
tions. The brass mount is fixed to a triangular massive 
metal plate supported by neoprene dampers. In our new 
setup, the neoprene damper positioned underneath the 
skimmer resides on a vertical translation stage and the 
remaining two dampers reside on two horizontal transla- 
tion stages. The translation stages are actuated by turn- 
ing knobs outside of the vacuum chamber which allows 
full adjustment of the nozzle position with respect to the 
beam axis even during operation (cf. Fig. [3]). 

In our present setup the maximum speed of rotation is 
limited to about 350 Hz (21000 rpm). The limiting factor 
so far is mechanical vibrations, which rapidly increase in 
amplitude at frequencies > 300 Hz. The origin of these 
vibrations is presumably a slight unbalance of the rotor 
as well as vibrations caused by the motor. These vibra- 
tions are found to reduce the lifetime of the ferrofluidic 
sealings of the feedthrough, which is specified to rotation 
frequencies < 166 Hz (10000 rpm). The motor as well 
as the CFK rotor, however, are able to withstand much 
higher frequencies up to 970 and 800 Hz, respectively. 



III. ANALYSIS OF TIME OF FLIGHT DATA 

Although the expansion out of the rotating nozzle is 
continuous, the beam passing through the skimmer and 
entering the detector chamber is chopped. Each time the 
nozzle comes close to the skimmer part of the ejected 
particles travel along trajectories that match the beam 
axis, determined by the positions of skimmer, quadrupole 
guides, and QMS detector. Using the vertical and hori- 
zontal translation stages the rotor position can be aligned 
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FIG. 3: Illustration of the dependence of measured signal 
(right side) on the alignment of the rotating nozzle setup with 
respect to the beam axis (dotted line), (a) Proper nozzle 
alignment; When the rotor is set to a perpendicular position 
with respect to the beam axis, nozzle and skimmer are in line 
with the beam axis up to the detector, (b) Nozzle shifted 
above the beam axis; Two pulses are produced each time the 
nozzle intersects the beam axis at two slightly tilted rotor 
positions, (c) Nozzle shifted below the beam axis; One pulse 
is measured with reduced amplitude and width. 



during rotation of the nozzle to optimize the detector 
signal. A peculiarity of the rotating-nozzle arrangement 
is the fact that the signal traces as a function of time 
are changed in an asymmetric manner as the rotor is 
shifted horizontally across to the beam axis, as illus- 
trated in Fig. [3J The signal traces are recorded using 
Ar gas expanding out of the nozzle at a nozzle pressure 
Po = 200mbar and at a frequency of rotation of 10 Hz. 
At this rotation frequency, the flight time of Ar atoms 
to the detector is short compared to the nozzle motion 
in front of the skimmer. Thus, the recorded signal traces 
mainly reflect the effective slit opening function given 
by the skimmer with respect to the moving nozzle. In 
Fig. [3J (a), the nozzle is properly aligned such that the 
signal trace features a slightly asymmetric trapezoidal 
shape. As the nozzle is moved across the beam axis by 
2 mm, the measured signal splits into two maxima, as 
shown in Fig. [3J (b). In this geometry, particle trajecto- 
ries matching the beam axis occur at two different angles 
of inclination of the rotor with respect to the beam axis. 



Rotor position (mm) 
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FIG. 4: Slit opening function measured at low rotation speeds 
of the nozzle (10 Hz) (solid line). The dashed line represents 
a model curve fitted to the data, consisting of the sum of 3 
Gaussian functions (light solid lines). 

Misalignment by 2 mm in the other direction leads to a 
reduction of the signal width and amplitude, as shown in 
Fig. [3J (c). The inclination signal envelopes in Fig. [3J (a) 
and (b) results from the decreasing solid angle defined 
by nozzle and skimmer aperture as the nozzle backs off 
away from the skimmer. 

In order to extract the exact values of the beam veloc- 
ity, temperature, and density from the measured signal 
traces, the effective slit opening function has to be char- 
acterized quantitatively and included into a model for 
the time-dependent density distribution of the molecular 
beam. The effective slit opening function for optimum 
horizontal nozzle alignment measured at low rotor fre- 
quency (10 Hz) is plotted in Fig. Q] (solid line). At such 
low rotor frequency longitudinal dispersion of the beam 
pulse during flight time is negligible. This is verified mea- 
suring the signal intensity at rest as a function of the ro- 
tor position using a micrometer driven manipulator. The 
resulting curve has identical shape compared to the time 
of flight curves at low rotor velocities. The dashed red 
line in Fig.[5]represents a fitted model function consisting 
of the sum of 3 Gaussian functions, 

3 (*-* f ) 2 

S2 (z) = ^A ie ~^r, (i) 

where amplitude factors Ai, center positions z i7 and 
widths Azi are determined by the fit to the experimental 
data. The density distribution in phase-space of a gas 
pulse produced by the rotating nozzle is modeled by the 
expression 

/o (r, v)= N f r (r ) f v (v), (2) 

which is normalized to the total number of molecules per 
bunch, N. 
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The spatial distribution f r is assumed to be anisotropic 
according to 



fr(r) = f±(x,y)f z (z). 
The distribution in transverse direction, 



(3) 



f±(x,y)=f ± e-—r 



(4) 

is normalized to fj_ — (7rAr 2 ) -1 . The distribution along 
the beam axis (z-axis) is obtained from the slit opening 
function (Eq. (JT)) ) by rescaling the z-coordinate according 
to z — > z/q, where q = \v s /v rot \ denotes the ratio between 
beam velocity v s = ^/k/(k — l)v and the velocity of the 
rotor tip, v ro t = 2irRf rot . Here, n — 5/3 for monoatomic 
gases, v = y / 2kBTo/m is the most probable velocity of 
the molecules inside the nozzle at a temperature To, and 
R = 19 cm stands for the rotor radius. Thus, 



/.(*) = f°/q 



3 



(5) 



where the amplitude factor f~ is given by f z = 
(y/Tr J2i=i Ai Az,;) _1 . The total number of molecules 
per pulse (integral over space of / r (r)) is a func- 
tion of speed of rotation, N — NAt, where N = 
A Po /(k B T ) v y/2 K /(K + 1)((k + l)/2) 1 /(i-«) [29] is the 
flux of molecules out of the nozzle orifice of area A — 
ir(d/2) 2 at pressure po and At = Az s n t /v ro t ~ 2cm/v rot . 
In other words, as the nozzle velocity v ro t falls below the 
beam velocity v s , the molecule bunches become longer 
than the spatial width of the slit opening function Az s m 
and the number of molecules per bunch increases accord- 
ingly. 

The velocity distribution is chosen according to the 
most frequently used ellipsoidal drifting Maxwellian 
model f v (v) = f vx (v x ,v y ) x f vz (v z ), where 



fv± — f v 



u± , and f vz = f vz e 



(6) 



with normalization constants f® = (wAv^) 1 and 
f° x = (i/nAvz)' 1 [H, The 1 widths of the veloc- 

ity distributions, Av± and Av z are related to the trans- 
verse and longitudinal beam temperatures T± and Tm by 
T|| ( j_ = mAwJ j_/(2fcs). While T\\ results from adiabatic 
cooling of the gas in the jet expansion, T± is an approxi- 
mate measure of the transverse beam divergence. Typical 
values of Tji attainable in free jet expansions range be- 
tween Tm ~ IK and T\\ ~ 10 K, while the transverse beam 



divergence corresponds to Tj_ ~ 300 K. During free ex- 
pansion, in the absence of collisions and external forces, 
the distribution function remains unchanged according 
to Liouville's theorem. The time evolution is given by 



f{r,v,t) = f (r- vt,v). 



(7) 



The spatial density as a function of time is then obtained 
by integrating over velocities, 



(8) 



n{r,t) = / d J vf(f,v,t). 




1,0 1,5 
Time of flight (ms) 



FIG. 5: Simulated density distributions on the beam axis as 
a function of time of flight at various speeds of rotation of the 
nozzle. Solid lines are calculated according to Eq. ((9} based 
on the measured slit-opening function, dashed lines reflect the 
longitudinal velocity distribution rescaled to flight times. 



The analytic solution of this expression is given by 



n(r,t) = Nnx(x,y)n z (z), 



where 



and 



3 

i=l 



Av 2 ± 



AvI 



Ar 2 



x +y 



Azt 



-1/2 



(9) 



(10) 



(11) 



Fig. [5] displays simulated Ar density distributions on 
the beam axis, n(0, 0, £, t), according to Eq. ((9]) for typical 
beam temperatures Tji = 5 K and T± = 323 K at various 
speeds of rotation (solid lines). Here, i = 0.38 m denotes 
the total flight distance from the nozzle to the QMS de- 
tector. The nozzle speed is v ro t — — 179 m/s (—150 Hz) 
in Fig. [5] (a), 24m/s (20 Hz) (b), 179m/s (150 Hz) (c), 
and 358 m/s (300 Hz) (d). Besides shifting linearly to 
longer flight times at higher rotational speeds, the distri- 
bution is first broadened at rotations around Hz due 
to the effect of the slit opening function. At higher 
speeds, the distribution first narrows down as the noz- 
zle moves faster past the skimmer (150 Hz), and then be- 
comes broader again as a consequence of longitudinal dis- 
persion, i. e. due to the fact that longer flight times allow 
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FIG. 6: (a) Measured time of flight density distributions of 
decelerated argon beams at various speeds of rotation of the 
nozzle, (b) Same data as in (a), where the horizontal axis is 
rescaled to velocities. 



the molecule bunches to longitudinally broaden due to 
the finite beam temperature Tji. The dashed lines repre- 
sent the neat velocity distributions horizontally rescaled 
to flight times according to t — £/v z . Thus, at nega- 
tive speeds of rotation (molecular beam is accelerated) 
the measured signal is predominantly determined by the 
shape of the velocity distribution, while at speeds around 
Hz the measured trace is strongly distorted by the slit- 
opening function. At higher rotation frequencies, the 
effect of the slit function diminishes and the two dis- 
tributions eventually converge. Therefore, at rotations 
/rot <; 200 Hz, Eq. (fTTj) may be replaced by the simpler 
expression n'(t) — A'/t 2 exp (— (s/t — vq) 2 /Av 2 ) in the 
procedure of fitting the measured data. 



IV. BEAM PARAMETERS 

Typical time of flight distributions measured with Ar 
at a nozzle pressure po = 400 mbar arc displayed in Fig. [6] 
(a) for various speeds of rotation. In these measurements, 
the Ar gas pressure inside the gas line was carefully ad- 
justed such that the effective pressure at the nozzle re- 



FIG. 7: Peak values of the time of flight argon density dis- 
tributions as a function of the rotation frequency. Symbols 
are the measured values, the solid line represents the result 
of fitting (Eq. (J9]l) to the data. 

mains constant. Due to the centrifugal force the nozzle 
pressure po rises sharply at increasing speed of rotation 
v ro t while holding the gas-line pressure po,gas line constant 
according to p = p , g as line exp (mv 2 ot / (2k B T )) [3, [2]. 
Signal broadening due to the convolution with the slit 
opening function around Hz as well as broadening due 
to longitudinal dispersion at high rotation frequencies are 
clearly visible (cf. inset of Fig. [6] (a)). The time of flight 
distributions shown in Fig. [5] (a) are horizontally rescaled 
to velocity distributions according to v = i/t and replot- 
ted in Fig. [6] (b). In this representation of the data the 
signal broadening due to dispersion at high rotation fre- 
quencies cancels out. However, broadening due to the 
slit opening function still remains and masks the real ve- 
locity distribution. The absolute densities are estimated 
by scaling the QMS detector signal to the output of the 
pressure gauge in the detector chamber. Besides, densi- 
ties can be estimated from the beam intensity measure- 
ment using the Pitot as discussed in Sec. [ill However, 
the density values obtained in this way turn out to be 
about a factor of 10 higher than the ones from the QMS 
measurement. This may be due to some feedback effect 
of the cold cathode ion gauge onto the vacuum inside the 
small volume of the Pitot tube, thus causing a nonlinear 
response of the vacuum gauge. 

Strikingly, the Ar signal amplitude rapidly drops by 
about a factor of 30 when increasing the speed of rota- 
tion from Hz to 300 Hz. In contrast, when the beam 
is accelerated by spinning the nozzle in the opposite di- 
rection up to 250 Hz, the signal only slightly increases 
by a factor of 2. The dependence of the measured peak 
density from the speed of rotation is shown in Fig. [7] as 
filled symbols. The solid line represents maximum val- 
ues of the model function, max(n), according to Eq. 
Here, the amplitude factor of the model function, N, is 
scaled by a factor of 0.5 to fit the experimental data. This 
strong amplitude modulation as a function of sense and 
speed of rotation is caused by both longitudinal as well 
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FIG. 8: Measured time of flight density distributions of kryp- 
ton at various speeds of rotation of the nozzle. 



as transverse beam dispersion. Since the rate of the sig- 
nal drop is determined by the beam temperatures Xj| and 
Tj_ , it is evident that beam temperatures have to be kept 
low by ensuring optimum free jet expansion conditions, 
i. e. high effective nozzle pressure po,nozz ^ 100 mbar at 
sufficiently low background pressure pb g < 10~ 3 mbar as 
well as a high quality nozzle aperture and skimmer (30| . 

Since the terminal beam velocity of jet beams is re- 
duced with increasing mass of the beam particles, decel- 
eration of beams of the heavier nobel gas atoms krypton 
(Kr) and xenon to small velocities is achieved at lower 
speeds of rotation. Beams of Kr, e.g., have terminal 
velocities around 400 m/s, such that deceleration down 
to velocities below 100 m/s are easily achieved at rota- 
tion frequencies > 250 Hz, as shown in Fig. [5J However, 
due to longitudinal dispersion, consecutive beam packets 
start to overlap at frequencies > 250 Hz ("wrap around 
effect" ) , which impedes quantitative analysis of the time 
of flight distributions. One way of passing around this 
effect could be the implementation of an additional me- 
chanical beam chopper to pick out individual beam pack- 
ets. 

In comparison with other methods of producing beams 
of slow molecules, our results are quite similar regard- 
ing beam densities at low velocities. Using electrostatic 
quadrupole filters, e.g., typical peak densities around 
10 9 cm -3 at velocities around 60 m/s are achieved [3l| . 
However, the rotating nozzle setup offers the possibility 
of tuning beam velocities in a wide range. Due to free 
jet expansion conditions, the internal degrees of freedom 
of the molecules are expected to be cooled to the Kelvin 
range, which can be reached using quadrupole filters only 
with substantial experimental efforts [32]. Besides, the 
rotating nozzle is suitable for any atomic of molecular 
substance available in the gas-phase at ambient temper- 
atures, in contrast to filtering or deceleration techniques, 
which rely on particularly suitable Stark of Zeeman ef- 
fects of molecules in external electrostatic or magnetic 
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FIG. 9: Argon and krypton beam parameters as a function of 
nozzle pressure. 



fields. 

In order to characterize the beam parameters and to 
test the analytic model, flight distributions of decelerated 
Ar and Kr beams are measured at various expansion pres- 
sures and speeds of rotation and are fitted by the model 
function, Eq. ©. Free fit parameters are the total num- 
ber of atoms per bunch, N, the reduced beam velocity vq, 
and the longitudinal beam temperature Tn . The result- 
ing values for Ar and Kr beams are depicted in Fig. [5] as 
symbols. Clearly, the results do not depend on the speed 
of rotation of the nozzle within the experimental uncer- 
tainty. As expected, the peak densities of Ar and Kr 
beams measured at rotation frequencies 150 and 100 Hz, 
respectively, evolve roughly proportionally to the nozzle 
pressure po (Fig. [5] (a)). The Kr densities fall behind the 
ones reached with Ar by about a factor of 3 due to the 
smaller flux of Kr out of the nozzle and due to the smaller 
longitudinal beam velocity of Kr. The fitted beam veloc- 
ities v s shown in Fig. [5] (b) are found to perfectly match 
the theoretical values for supersonic beams, v s , provided 
the nozzle pressure exceeds about 100 mbar. The latter 



are given by v s 



2k ksTo 
K—l m 



1 



( P^] 



where Tn 



and p are the nozzle temperature and pressure, pb g is the 
background gas pressure and k = cp/cy = 1 + 2/ / is the 
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Nozzle pressure (bar) 

FIG. 10: Attenuation of the krypton peak beam intensity as 
a function of nozzle pressure for different speeds of rotation 
of the nozzle. 



adiabatic exponent for a gas of particles with / degrees 
of freedom [29|]. For mono-atomic gas (k = 5/3) and low 
background pressure pb g <C po we obtain the simpler for- 
mula v s — y/5 kBTo/m, which yields v s — 579 m/s for 
Ar and v s = 399 m/s for Kr at T = 323 K. At lower 
nozzle pressure pg , the fitted values of Vq slightly deviate 
from the model curve, indicating expansion conditions at 
the transition from the effusive to the free jet expansion 
regime. 

Fit values of the longitudinal temperature Tji of Ar 
beams are shown in Fig. [5] (c) as symbols. The solid 
lines represent model curves according to two differ- 
ent collision models. Assuming hard-spheres (HS) col- 
lisions, the terminal speed ratio S — v s / Av is given by 
S = 0.289(poda /Tq) oa , while atom-atom interactions ac- 
cording to the Lennard-Jones (LJ) potential yield S = 

0.397(p de 1/3 ^/T 4/3 ) - 53 |2p. Here, a = 46.3 [50.3] A 2 
represents the HS-cross section, e = 141.6 [50.3] Kxfce 
and r m — 3.76 [4.01] A stand for the potential well depth 
and the equilibrium distance of the LJ-potential for argon 
[Kr], respectively. In comparison with the experimental 
data, the model assuming interactions according to the 
LJ-potential appears to yield more precise values than 
the one assuming HS-collisions. However, the LJ-model 
systematically slightly under-estimates longitudinal tem- 
peratures, in contrast to the HS-model, which gives too 
large values. The measured temperature continuously 
decreases as a nozzle pressure up to about 1 bar is ap- 
plied, reaching values T» < 1.5 K with both Ar and Kr. 
Evidently, high nozzle pressures are desirable if low tem- 
peratures are to be achieved. 



V. DETRIMENTAL EFFECTS 

In order to work at high nozzle pressures while keep- 
ing the gas load low it may be favorable to use smaller 
apertures, which goes at the expense of beam intensity, 



though. The maximum applicable gas pressure, how- 
ever, is limited by the pumping speed available in the 
rotor chamber. In our setup, the rotor-chamber pres- 
sure increases roughly linearly from 1.9 x 10~ 5 mbar to 
2.2 x 10 _4 mbar, as the Ar pressure at the nozzle rises 
from lOOmbar up to 1000 mbar. At these vacuum con- 
ditions, collisions of beam molecules with the ones from 
the residual gas may seriously deplete the beam. This ef- 
fect is studied by measuring the peak beam density of Kr 
beams as a function of nozzle pressure up to po < 2 bar at 
different rotation frequencies, / = 50 Hz (vq = 340 m/s), 
/ = 100 Hz (v = 280 m/s), / = 150 Hz (v = 220 m/s). 

As shown in Fig. [TU1 the Kr beam density is affected 
by rest gas collisions as the nozzle pressure exceeds about 
600 mbar. Nozzle pressures higher than about 900 mbar 
even lead to the destruction of the beam and beam densi- 
ties drop down again, irrespective of the rotor frequency. 
Only at high frequencies leading to beam velocities be- 
low about 100 m/s the position of the maximum of the 
dependence n(po) shifts down to lower values of po due to 
the fact that slow molecules spend more time inside the 
region of high rest-gas pressure and are therefore more 
sensitive to rest-gas collisions. The solid lines in Fig. [10] 
are fit curves according to the model of the peak den- 
sity (Eq. ©), taking into account a density reduction 
due to collisions, n re d = nexp (— arib g v re i s/vq). Here, 
the relative velocity between colliding atoms is taken as 

v re i — tj Vq + v% g , s ~ 2 cm is the average flight dis- 
tance from nozzle to skimmer, and nb g = Pbg/ksT and 
Vf,g = y^Sk^To/ (wm) are the residual gas density and 
mean thermal velocity, respectively. The Kr-Kr collision 
cross section a is varied as a free fit parameter, yield- 
ing an average value a = 130 ± 30 A . This value is m 
reasonable agreement with the Landau-Lifshitz approxi- 
mation for the maximum scattering cross section, <j max = 

8.083 (C^/(hv re i)) V5 = 390A 2 , with = 77.4 eVA 6 
for Kr [H[. 

High nozzle pressures, which lead to increased rest-gas 
pressure in the nozzle chamber, cause another unwanted 
effect to occur - the formation of a continuous effusive 
beam. The latter adds up to the pulsed nozzle beam 
and may have comparable beam intensities as the nozzle 
beam is decelerated to low speeds and low densities. We 
therefore compare the densities of nozzle beams and effu- 
sive beams by recording the detector signal with the noz- 
zle standing properly aligned in front of the skimmer, on 
the one hand, and with the rotor being tilted completely 
out of the beam axis, on the other. The resulting density 
values are plotted in Fig.QT] While the measurement us- 
ing the Pitot tube yields higher nozzle beam densities by 
a factor of 10, equal densities are found for the effusive 
beam. Thus, the nozzle beam is higher in density by 3-4 
orders of magnitude than the effusive beam. 

For decelerated beams of Ar and Kr, however, equal 
densities of nozzle and effusive beams are reached at 
beam velocities vo < 60 and vq < 40 m/s, respectively. 
Provided the expansion is operated in the supersonic 
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FIG. 11: Comparison of beam densities of the continuous noz- 
zle beam (rotor at rest) and the effusive beam produced by 
the residual gas in the rotor chamber. Circles depict the data 
measured using the quadrupole mass spectrometer scaled to 
densities obtained from the detector chamber pressure mea- 
sured using a cold cathode pressure gauge. Squares represent 
the density values obtained from the Pitot tube measurement. 

regime, the ratio of nozzle and effusive beam densities 
is independent of the nozzle pressure po- This ratio can 
be improved either by increasing the pumping speed in 
the nozzle chamber, which is characteristic for such an 
arrangement based on a continuous gas expansion. The 
effusive background can be further suppressed by imple- 
menting a mechanical velocity selector synchronized with 
the rotating nozzle. Therefore, in the near future a chop- 
per wheel will be placed behind the skimmer. 



VI. ELECTROSTATIC GUIDING 

The issue of transverse divergence of decelerated beams 
can be relieved for beams of polar molecules by im- 
plementing additional guiding elements using inhomoge- 
neous electrostatic fields. This technique has been used 
in the past for focussing and even state selecting molecu- 
lar beams for reactions dynamics experiments [3J] . More 
recently, the technique has been used for guiding and 
trapping of cold polar molecules produced by Stark de- 
celeration or Stark filtering 0,[Ia[H, Ell- 
in our setup we have implemented an electrostatic 
quadrupole guide between the skimmer and the ionizer 
of the QMS. It consists of four 259 mm long stainless 
steel rods 2 mm in diameter with a gap of 2 mm between 
two diagonally opposing rods. In a first experiment, slow 
beams of fluoroform (CHF3) molecules have been pro- 
duced by expanding CHF 3 either as a neat gas or seeded 
in Kr out of the rotating nozzle. CHF3 is a polar molecule 
with a permanent dipole moment of 1.65Debye. Since 
its mass (mcHF3 ~ 70amu) lies close to the mass of Kr 
(rriKr ~ 80 amu) it can be efficiently cooled and deceler- 



FIG. 12: Typical CHF3 density distributions as a function of 
flight time for two different speeds of rotation of the nozzle. 
The data recorded at 250 Hz are scaled up by a factor of 
400. Signals shown as back lines are recorded with no voltage 
applied to the guide electrodes, light red lines are obtained at 
3.5 kV electrode voltage. 



ated by coexpansion with Kr. 

When applying a high voltage U = 3.5 kV to a pair 
of opposing electrodes the QMS detector signal clearly 
rises, as depicted in Fig. [T2l In these measurements, 
CHF3 is expanded at a nozzle pressure po = 200 mbar at 
rotor frequencies 75 and 250 Hz, leading to most probable 
beam velocities v w — 450 and v w — 165 m/s, respectively. 
Clearly, the enhancement of the QMS signal at an elec- 
trode voltage U — 3.5 kV in proportion to the signal with 
no high voltage is higher at low beam velocities, reaching 
roughly a factor of 5 for v w = 165 m/s. Great care was 
taken to suppress any direct influence of the electrode 
voltage on the ionization efficiency of the QMS ionizer. 

Similar experiments have also been performed with 
deuterated ammonia (ND 3 ). In addition to increasing 
the beam intensity, the quadrupole fields act as quan- 
tum state selectors which filter out those states that are 
attracted towards the beam axis where the electric field 
amplitude is minimal (low-field seeking states). Provided 
the rotational temperature of the ND3 molecules is com- 
parable to the translational one, i. e. in the range of a 
few Kelvin, we may expect to produce ND3 molecules in 
the state J = 1, M = K = —1 with very high purity [35| . 

Interestingly, the dependence of the enhancement of 
peak density on the electrode voltage U features a non- 
monotonic increase, as shown in Fig. 1131 A first local 
signal maximum is observed at electrode voltages be- 
tween 200 and 500 V, depending on beam velocity. We 
interpret this behavior in terms of two competing mech- 
anisms being active: Guiding of molecules in low-field 
seeking states enhances the detector signal whereas devi- 
ation of high-field seeking molecules away from the beam 
axis reduces the measured beam density. At low voltages 
U ~ 500 V, the guiding effect outreaches the deviation 
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The dependence of the enhancement factor on the 
beam velocity is illustrated in Fig. [TJ] In this mea- 
surement, the electrode voltage is set to the maximum 
value U = 3.5 kV and the beam velocity is tuned from 
v w = 500 down to v w = 250 m/s by changing the fre- 
quency of rotation of the nozzle. Further deceleration 
down to v w « 150 m/s is obtained by seeded expansion 
of CHF3 in Kr at a mixing ratio 1:2. The enhancement 
factor continuously rises as the beam velocity is reduced, 
which highlights the potential of combining decelerated 
nozzle beams with electrostatic guiding fields. 



FIG. 13: Relative efficiency of guiding CHF3 molecules at 
various beam velocities as a function of electrode voltage. 



VII. CONCLUSION AND OUTLOOK 
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FIG. 14: Enhancement of the CHF3 peak density as a func- 
tion of beam velocity for fixed electrode voltage U = 3.5 kV. 



efficiency, leading to a local density maximum. As the 
electrode voltage rises, all high-field seeking molecules are 
expelled out of the beam and only low-field seekers are 
guided to the detector. This guiding effect partly com- 
pensates the transverse blowing up of the beam, which is 
most relevant for low longitudinal beam velocities. Con- 
sequently, the relative guiding efficiency increases up to 
a factor of 5 as the beam velocity is reduced down to 
v w = 165 m/s. The fact that the same qualitative behav- 
ior is observed with beams of ND3 supports our interpre- 
tation, which does not consider specific Stark shifts of the 
molecules. In future efforts the transmission of our setup 
will be further studied both experimentally and using 
molecular trajectory simulations. Besides, guiding of de- 
celerated molecules around bent electrodes will be imple- 
mented to provide state-purified beams of cold molecules. 
This has the additional advantage of efficiently suppress- 
ing the hot effusive background. 



In conclusion, we have presented a versatile appara- 
tus that produces cold beams of accelerated or deceler- 
ated molecules based on a rapidly rotating nozzle. Var- 
ious technical improvements with respect to the original 
demonstration by Gupta and Herschbach p], Q are intro- 
duced. In particular, gas injection into the rotor is now 
realized using a fluid metal-sealed rotary feedthrough 
which eliminates gas leakage into the vacuum chamber. 
Using this setup, beam velocities well below 100 m/s and 
longitudinal beam temperatures down to about 1 K are 
achieved. A detailed analytic model to simulate the mea- 
sured density distributions as functions of flight time is 
presented. The fundamental drawback of this technique, 
the sharp drop of beam intensity at slow beam veloc- 
ities, is relieved by combining the setup with electro- 
static guiding elements, provided the molecules feature 
a suitable Stark effect. Moreover, Stark guiding goes 
along with internal-state selectivity. Thus, a significant 
increase of beam intensity as well as an expected high 
state-purity is achieved with beams of fluoroform and 
ammonia molecules. 

Thanks to the technical improvements the rotating 
nozzle setup has evolved into a simple and reliable source 
of cold molecular beams with tunable velocity. In partic- 
ular, slow beams of molecules that are not amenable to 
efficient filtering or deceleration by means of the Stark 
interaction with external electric fields can now be pro- 
duced with high intensity. In the near future, we will uti- 
lize this beam source for studying reactive collisions be- 
tween slow molecules and cold atoms produced by other 
means, e. g. using a magneto-optical trap. At collision 
energies down to < 1 meV, interesting new quantum ef- 
fects may be expected in the reaction dynamics, e. g. 
strong modulations of the scattering cross-section in re- 
actions of the type Li+HF^LiF+H 27]. Furthermore, 
intense cold molecular beams with tunable velocity and 
internal-state purity are of current interest for surface 
scattering experiments 28] . 
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